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A simple and accurate approximation for the Q stability 
parameter in multi-component and realistically thick discs 
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ABSTRACT 

In this paper, we propose a Q stability parameter that is more reahstic than those 
commonly used, and is easy to evaluate [see Eq. (19)]. Using our Qn parameter, you 
can take into account several stellar and/or gaseous components as well as the stabi- 
lizing effect of disc thickness, you can predict which component dominates the local 
stability level, and you can do all that simply and accurately. To illustrate the strength 
of Qat, we analyse the stability of a large sample of spirals from The H i Nearby Galaxy 
Survey (THINGS), treating stars. Hi and H2 as three distinct components. Our anal- 
ysis shows that H2 plays a significant role in disc (in) stability even at distances as 
large as half the optical radius. This is an important aspect of the problem, which was 
missed by previous (two-component) analyses of THINGS spirals. We also show that 
H I plays a negligible role up to the edge of the optical disc; and that the stability level 
of THINGS spirals is, on average, remarkably flat and well above unity. 

Key words: instabilities - stars: kinematics and dynamics - ISM: kinematics and 
dynamics - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: star forma- 
tion. 



1 INTRODUCTION 

Today, several decades after the pioneering work of Lin 
& Shu (1966) and the seminal papers by Jog & Solomon 
(1984a, b), it is widely accepted that stars and cold inter- 
stellar gas have an important interplay in the gravitational 
instability of galactic discs. The gravitational coupling be- 
tween stars and gas does not alter the form of the local 
axisymmetric stability criterion, Q > 1 (Toomre 1964), but 
makes the Q stability parameter dependent on the radial 
velocity dispersions and surface mass densities of the two 
components (Bertin & Romeo 1988; Elmegreen 1995; Jog 
1996; Rafikov 2001; Shen & Lou 2003). The value of Q is 
also affected by other factors, such as the vertical structure 
of the disc (Shu 1968; Romeo 1990, 1992, 1994; Elmegreen 
2011; Romeo & Wiegert 2011), gas turbulence (Hoffmann 
& Romeo 2012; Shadmehri & Khajenabi 2012), and gas 
dissipation (Elmegreen 2011). Comprehensive analyses have 
shown that the two-component Q parameter has a large im- 
pact on the dynamics and evolution of spiral structure in 
galaxies (see Bertin & Lin 1996), and is also a useful di- 
agnostic for exploring the link between disc instability and 
star formation (Leroy et al. 2008). 

Romeo & Wiegert (2011) introduced a simple and ac- 
curate approximation for the two-component Q parameter. 
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which takes into account the stabilizing effect of disc thick- 
ness and predicts whether the local stability level is dom- 
inated by stars or gas. The Romeo- Wiegert approximation 
has been used for investigating the evolution of gravitation- 
ally unstable discs (Cacciato et al. 2012; Forbes et al. 2012), 
the spiral structure of NGC 5247 (Khoperskov et al. 2012), 
the dynamical link between dark matter and Hi in nearby 
galaxies (Meurer et al. 2013), as well as the link between disc 
stability and the relative distributions of stars, gas and star 
formation (Zheng et al. 2013). Forbes et al. (2012) concluded 
that the Romeo- Wiegert approximation is much faster to use 
than the Q stability parameter of Rafikov (2001): it speeds 
up their disc-evolution code by as much as one or two orders 
of magnitude! This is simply because such an approximation 
estimates Q analytically, without the need to minimize the 
dispersion relation over all wavenumbers, as is usually done. 

A fundamental problem that must be faced when 
analysing the stability of galactic discs is how to represent 
their complex structure using only two components. The re- 
sults of such analyses are indeed very sensitive to the choice 
of the gaseous ID velocity dispersion, a^: the colder the gas, 
the stronger its impact on the stability of the disc (e.g.. 
Jog & Solomon 1984a; Bertin & Romeo 1988). Choosing 
cTg « 6 kms~^ will represent molecular gas well (Wilson et 
al. 2011), but will overestimate the contribution of atomic 
gas and make the disc more unstable than it actually is. 
Vice versa, choosing Ug « 11 kms~^ will represent Hi well 
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(Leroy et al. 2008), but will underestimate the contribution 
of H2 and stabilize the disc artificially. Intermediate values 
or more elaborate choices of (jg will still not solve the prob- 
lem. This motivates the use of a proper multi-component Q 
parameter. 

One of the first papers that discussed the gravitational 
instability of multi-component discs dates back to Moro- 
zov (1981). This author derived a dispersion relation that is 
valid for infinitesimally thin discs made of gas and TV* stel- 
lar components. He also calculated the stability criterion for 
A'* = 2. Rafikov (2001) derived a stability criterion that is 
valid for any A^* and can easily be expressed in the usual 
form Q > 10 

In this paper, we introduce a new Q stability parame- 
ter, which has the same strong advantages as the Romeo- 
Wiegert approximation, and which is applicable to fully 
multi-component and realistically thick discs [see Eq. (19)]. 
We also show how to use our Qn parameter for analysing 
the stability of galactic discs, and illustrate the strength of 
a multi-component analysis. We do so using a large sample 
of spirals from The Hi Nearby Galaxy Survey (THINGS), 
and treating stars. Hi and H2 as three distinct components. 

The rest of the paper is organized as follows. In Sect. 

2, we review the basic case of two-component discs and fur- 
ther motivate the need for a multi-component analysis (Sect. 
2.1), we present our Qm parameter (Sects 2.2-2.4), and we 
analyse the stability of THINGS spirals (Sect. 2.5). In Sect. 

3, we discuss the weaknesses of Qn, which are common to 
all Q parameters and stability criteria quoted here. In Sect. 

4, we draw the conclusions. 



2 STAR-GAS INSTABILITIES AND THE Q 
DIAGNOSTIC 

2.1 Two components ... or more? More! 

Let us first discuss the case of two-component and infinites- 
imally thin discs, which is fundamental to a proper under- 
standing of Sects 2.2-2.5. It is well known that the stability 
properties of such discs are determined by five basic quanti- 
ties: the epicyclic frequency, k, the stellar and gaseous sur- 
face densities, and Eg, the stellar radial velocity disper- 
sion, (Jj;*, and the gaseous ID velocity dispersion, criDg (e.g., 
Lin & Shu 1966; Rafikov 2001). The last two quantities re- 
flect an important dynamical difference between stars and 
cold interstellar gas. The stellar component is collisionless so 
its velocity dispersion is anisotropic, while gas is coUisional 
and has an isotropic velocity dispersion (see, e.g., Binney & 
Tremaine 2008) . Hereafter we will simplify the notation and 
denote the relevant velocity dispersions of the two compo- 
nents with (J* and erg. 

Lin & Shu (1966) and Rafikov (2001) took those facts 
into account by treating stars as a kinetic component and 
gas as a fiuid (see also Shu 1968). The resulting local ax- 
isymmetric stability criteria are equivalent because they are 
based on the same dispersion relation, ct)'^(fc), and because 
they are derived by imposing that uj'^ (k) > for all k. How- 
ever, Rafikov's stability criterion is simpler and more used 



^ Another (unpublished) stability analysis of A'^-component discs 
was made by Romeo (1985), pp. 140-145 and 215-216. 
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Figure 1. Accuracy of the two- fluid stability criterion. The 
curves shown are the contour lines of the relative error (Qnf — 
QRk)/QRki where Qj^f and Qjik a^re the fluid-fluid and kinetic- 
fluid Q parameters of Rafikov (2001). In addition, and Qg are 
the stellar and gaseous Toomre parameters, tr* and Cg are the 
relevant velocity dispersions of the two components. 



(e.g., Dalcanton et al. 2004; Li et al. 2005, 2006; Kim & 
Ostriker 2007; Yang et al. 2007; Yim et al. 2011). Such a 
criterion can be written as Q^k > 1, where 



— — = max{SCRk(Ar)} , 
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In these equations, Qbm is Rafikov's kinetic-fluid Q param- 
eter, SCRk(i^) is the related stability curve, K is the radial 
wavenumber of the perturbation expressed in dimensionless 
form, Qi, and Qg are the stellar and gaseous Toomre pa- 
rameters, and Jo denotes the modifled Bessel function of the 
first kind and order zero. 

Berlin & Romeo (1988), Elmegreen (1995), Jog (1996), 
and again Rafikov (2001) adopted a less rigorous, but more 
straightforward approach: they treated the stellar compo- 
nent as a fluid, with sound speed equal to cr* . The resulting 
two-fluid stability criteria are equivalent, apart from the dif- 
ferent parametrizations used, because they are based on the 
same dispersion relation (Jog & Solomon 1984a). However, 
even in this case, Raflkov's criterion is simpler and is be- 
coming more and more widely used (e.g., Leroy et al. 2008; 
Robertson & Kravtsov 2008; Dekel et al. 2009; Mastropietro 
et al. 2009; Ceverino et al. 2010; WestfaU et al. 2011; Elson 
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Figure 2. Accuracy of the Romeo- Wiegert approximation. The 
curves shown are the contour lines of the relative error (Qrw ~ 
QRk)/QRki where Qrw is the two-component Q parameter of 
Romeo & Wiegert (2011), and QRk is the kinetic-fluid Q param- 
eter of Rafikov (2001). The rest of the notation is the same as in 
Fig. 1. 



et al. 2012; Watson at al. 2012; Williamson & Thacker 2012). 
Such a criterion can be written as Qnf > 1, where Rafikov's 
fluid-fluid Q parameter can be computed by maximizing the 
related stability curve over all radial wavenumbers: 



1 



:max{SCRf(i^)}. 
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(6) 
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What is the accuracy of the two-fluid stability criterion? 
Bertin & Romeo (1988) compared the fluid-fluid and kinetic- 
fluid marginal stability curves in a representative set of 
cases, and showed that the differences are small (see their 
flg. 2). Rafikov (2001) carried out a more detailed analysis 
and got similar results (see his fig. 3). However, he did not 
evaluate the relative error (QRf — QRk)/QRk, which depends 
on s and 



(8) 



and which is a useful piece of information for assessing the 
accuracy of Qrj. We do this in Fig. 1. Our contour map 
shows that such error is less than 7% and has a root-mean- 
square value of 4%. This means that Q^s and Qnk are practi- 
cally equivalent, since the parameters s and q are themselves 
subject to observational uncertainties. 

Romeo & Wiegert (2011) showed that the two-fluid Q 
parameter can be accurately estimated without the usual 
maximization (or minimization) procedure: 
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Although very recent, such an approximation is already fre- 
quently used (e.g., Cacciato et al. 2012; Forbes et al. 2012; 
Hoffmann & Romeo 2012; Khoperskov et al. 2012; Meurer 
et al. 2013; Zheng et al. 2013). Fig. 2 shows that Qrw is 
remarkably accurate even with respect to the kinetic-fluid 
Q parameter. In fact, the relative error (Qrw — QRk)/QRk 
is below 10% and has a root-mean-square value of 5%. Thus 
Qrw is a faster and physically equivalent alternative to QRf 
or QRk. 

To understand the weaknesses of a two-component anal- 
ysis, let us see how spiral galaxies populate the parameter 
plane of star-gas instabilities. We use a sample of twelve 
nearby star- forming spirals from THINGS, previously anal- 
ysed by Leroy et al. (2008) and Romeo & Wiegert (2011): 
NGC 628, 2841, 3184, 3198, 3351, 3521, 3627, 4736, 5055, 
5194, 6946 and 7331. For each galaxy of this sample, we 
compute the radial profiles s = s{R) and q — q{R), and 
hence the track left by the galaxy in the {s,q) plane. The 
result for the whole sample is shown in the left panel of Fig. 
3. Note that 20% of the data fall within the shaded part 
of the (s, q) plane, which represents the condition for star- 
gas decoupling. In this region, SCRk(-ft^) has two maxima: 
one at small K, where the response of the stellar compo- 
nent peaks; and the other at large K, where gas dominates. 
In the 'stellar phase', the maximum at small K is higher 
than the other one, and therefore it controls the onset of 
disc instability. Vice versa, in the 'gaseous phase', it is the 
maximum at large K that determines QRk. The two-fiuid 
counterpart of this region (thin dashed lines) is also popu- 
lated by 20% of the data (Romeo & Wiegert 2011). In the 
rest of the parameter plane, the dynamical responses of the 
two components are strongly coupled and peak at a single 
wavelength. 

The analysis of THINGS spirals carried out above treats 
the interstellar medium (ISM) as a single component with 



Ehi + Sh2 and 



11 kms"^ (Leroy et al. 2008). 



What are the limitations of this approach? How do H I and 
H2 contribute to star-gas instabilities? To answer these ques- 
tions, we consider Hi and H2 separately, and choose obser- 
vationally motivated values of the ID velocity dispersion: 
(Thi = 11 kms~^ (Leroy et al. 2008), and aH2 = 6 kms~^ 
(Wilson et al. 2011). We then compute the {s,q) tracks for 
each case (stars plus Hi or H2), and show the results in the 
right panel of Fig. 3. Note that Hi and H2 populate the pa- 
rameter plane differently. In particular, none of the H I data 
falls within the two-phase region, while H2 populates such a 
region in 60% of the cases. This means that H I and H2 have 
distinct stability properties, and a fundamentally different 
dynamical coupling with stars. Treating the ISM as a single 
component underestimates the role that H2 plays in star-gas 
instabilities, and overestimates the contribution of Hi. This 
is why a multi-component analysis is needed! 
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Figure 3. Why a multi-component analysis is needed: the parameter plane of star-gas instabilities populated by THINGS spirals. The 
cases illustrated in the two panels are stars plus gas (left), and stars plus Hi or H2 (right). The galaxy data are from Leroy et al. (2008), 
Qi, and Qg are the stellar and gaseous Toomre parameters, cr* and Cg are the relevant velocity dispersions of the two components. The 
shaded part of the (s, g) plane represents the condition for star-gas decoupling, with the stellar component treated as coUisionless (thick 
solid lines) or coUisional (thin dashed lines). This is the 'two-phase region' discussed in the text. The boundaries of this region and the 
transition line intersect at (s, q) ~ (0.21, 0.88) in the coUisionless case, and at (s, q) ~ (0.17, 1) in the coUisional case. 



2.2 Approximating Q in the thin-disc hmit 

The only multi-component stability diagnostics that have 
been available so far are the local axisymmetric stability 
criteria of Morozov (1981) and Rafikov (2001), which are 
valid for infinitesimally thin discs made of gas and A'^* stel- 
lar components. For A''* = 2, the case considered by Mo- 
rozov (1981), such criteria are equivalent because they are 
based on similar approximations. However, besides being 
more general, Rafikov's criterion is simpler and can be writ- 
ten as Qr.jv > 1, where 



max{SCR,]v(fc)} 



SCR,iv(fc) = 7^ 



N 
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Let us illustrate how to find a simple, accurate, fast 
and more general Q diagnostic. In Sect. 2.1, we have shown 
that stars can be accurately treated as a fluid when evalu- 
ating Q. So we can safely replace the kinetic terms in Eq. 
(12) with their fluid counterparts. We then face the heart 
of the problem: how to estimate the least stable wavenum- 
ber, fcmax, without the usual maximization procedure. Con- 
sider the two-component case first, and compare the Romeo- 
Wiegert approximation [Eqs (9) and (10)] with the two-fiuid 
stability parameter [Eqs (6) and (7)]. One can easily infer 
that A'max fulfils the following conditions: 



(i) If > Qg, then 2K^^^/{1 + Kl,^) ~ W and 
2Km^^s/{l + K^^^s'^) ~ 1. As = 2s/{l + s'^), this implies 

that A'max ~ 1/s, i.e. femax ~ K/fTg . 

(ii) If Qg > Q*, then 2A:„,ax/(l + i^max) ~ 1 and 

2i^maxS/(l + i^maxS^) ~ W. This implies that A"max ~ 1, 
i.e. fcmax ^ ti/o'-k- 

Conditions (i) and (ii) have not been pointed out in previous 
analyses. They simply mean that 1/femax is approximately 
the typical epicycle size of the less stable component. This 
approximation is not accurate when s « 0.2 and q ~ 1, since 
for such values there is a transition between three stability 
regimes and A"max has a jump across one of the interfaces 
(see discussion of Fig. 3). Note, however, that estimating 
^i^max as above produces a very accurate estimate of Qm, 
namely the Romeo- Wiegert approximation. This is because 
Qjii is continuous across the q — 1 line, and because the 
error that affects the estimate of Qb.! is of second order with 
respect to that of ifmax: AQnf oc (Aifmax)^- The first-order 
term is obviously zero, since the first derivative of SCRf(if) 
vanishes for K = A'maxQ This flow of arguments suggests 
that we can estimate the A'^-component Q parameter as in 
the two-fluid case, i.e. by approximating 1/fcmax with the 
typical epicycle size of the least stable component (fcmax ~ 

ft/cTm): 



_L - V — 

Qat ^ ^ O. 



(15) 



^ This is actually the idea behind all minimization (or maximiza- 
tion) problems, and the reason why their solutions are robust. 
An instructive example is the 'optimal' Wiener filter used in sig- 
nal/image processing (see, e.g., Press et al. 1992). 
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Table 1. Accuracy of our approximation in the ten-component 
cases analysed by Rafikov (2001). 



Case 


Qr,n 


Qn 


(Qn — Qr,n)/Qr,n 


1 


1.10 


1.12 


1.8% 


2 


1.00 


1.04 


3.7% 


3 


1.00 


1.04 


4.5% 


4 


1.00 


1.01 


0.8% 


4+1 


1.00 


1.01 


1.2% 



Qb.,n is the A'^-component Q parameter of Rafikov (2001), Qiv is 
our TV-component Q parameter, and (Qjv — Qr,jv)/Qr,jv is the 
relative error of our approximation. 



Wi 



2a„ 



(16) 



where the index m denotes the component with smallest Q: 
Qm ~ minjQi}. This is the component that dominates the 
local stability level {Wm = 1). All other components have 
less weight; the more ai differs from am, the smaller the 
weight factor Wi. 

What is the accuracy of Qjv? In Appendix A, we show 
how estimation uncertainties propagate from Wi to Qjv , and 
derive an upper bound for the resulting root-mean-square 
error: 



Qjv 



< 0.03 VN 



(17) 



Eq. (17) tells us that the accuracy of our approximation de- 
teriorates slowly as we consider more and more components. 
In the case of greatest interest, stars plus Hi plus H2, the 
relative error is on average less than 6%, i.e. almost as low 
as in the two-component case (see Fig. 2). 

To demonstrate the accuracy of our approximation, we 
put it to a stringent test: the 10-component model of the So- 
lar neighbourhood analysed by Rafikov (2001). That model 
includes cold interstellar gas, giants, stars in six luminos- 
ity ranges, white and brown dwarfs. Rafikov also analysed 
the effect of varying the most uncertain model parameters. 
He motivated and discussed five cases, which we summa- 
rize in Appendix B. For each case, we compare Qn with 
Rafikov's Qr.jv parameter and compute the relative error 
{Qn — Qr,n)/Qr,n ■ Table 1 shows that our approximation 
is remarkably accurate. In all cases, the relative error is less 
than 5%, which means twice as low as the upper bound pre- 
dicted by Eq. (17). 



2.3 Adding the effect of disc thickness 

Our approximation is not yet complete. It does not include 
the stabilizing effect of disc thickness, which is important 
and should be taken into account when analysing the stabil- 
ity of galactic discs (Romeo & Wiegert 2011). In this section, 
we generalize the Romeo- Wiegert approach and provide a 
simple recipe for adding such an effect. 

From the thin-disc limit, we have learned that the local 
stability level is dominated by the component with smallest 
Q [see Eqs (15) and (16)]. The contributions of the other 
components are weakened by the Wi factors, which are dif- 
ferent and small if the components are dynamically distinct. 
In this case, we can estimate the effect of thickness reason- 



ably well by considering each component separately. Romeo 
(1994) analysed this case in detail. The effect of thickness 
is to increase the stability parameter of each component by 
a factor T, which depends on the ratio of vertical to radial 
velocity dispersion: 



l + O.ef — V for 0<a^/afl <0.5, 



0.8 + 0.7 f—") for 0.5<cr^/afl<l. 



(18) 



Eq. (18) can be inferred from fig. 3 (top) of Romeo (1994). 
The range 0<az/aR<0.5 is characteristic of the old stel- 
lar disc in Sc-Sd galaxies (Gerssen & Shapiro Griffin 
2012), while 0.5 < o-z/(t_r < 1 is the usual range of velocity 
anisotropy (typical of the old stellar disc in Sa-Sbc galaxies, 
of young stars and the ISM). To approximate Q in this more 
general context, use then Eq. (15) with Qi replaced by TiQi: 



1 



JV 

1=1 



(19) 



where Qjv is our Q stability parameter for multi-component 
and realistically thick discs, Qi = ftai/vrGEi is the Toomre 
parameter of component i, Ti is given by Eq. (18), and Wi is 
given by Eq. (16). NOTE that the index m now denotes the 
component with smallest TQ: TmQm ~ min{TiQi}. This is 
the component that dominates the local stability level. The 
contributions of the other components are still suppressed 
by the Wi factors. 

2.4 What about the effect of ISM turbulence? 

Turbulence plays a fundamental role in the dynamics and 
structure of cold interstellar gas (see, e.g., Elmegreen & 
Scalo 2004; McKee & Ostriker 2007; Agertz et al. 2009). The 
most basic aspect of interstellar turbulence is the presence of 
supersonic motions. These are usually taken into account by 
identifying erg with the typical ID velocity dispersion of the 
medium, rather than with its thermal sound speed. Another 
important aspect of interstellar turbulence is the existence 
of scaling relations between Eg, cjg and the size of the region 
over which such quantities are measured {£). Observations 
show that Ehi ~ £^^^ and am ~ ^^''^ up to scales of 1-10 
kpc, whereas Eh2 ~ constant and (jh2 ~ £^^^ up to scales 
of about 100 pc (see, e.g., Elmegreen & Scalo 2004; McKee 
& Ostriker 2007; Romeo et al. 2010). 

Motivated by the large observational uncertainties of 
Eg(£) and ag{£), and having in mind near- future applications 
to high-redshift galaxies, Romeo et al. (2010) considered 
more general scaling relations. Eg oc and ag cx £'', and 
explored the effect of turbulence on the gravitational insta- 
bility of gas discs. They showed that turbulence excites a rich 
variety of stability regimes, several of which have no classi- 
cal counterpart. See in particular the ^stability map of tur- 
bulence' (fig. 1 of Romeo et al. 2010), which illustrates such 
stability regimes and populates them with observations, sim- 
ulations and models of ISM turbulence. Hoffmann & Romeo 
(2012) extended this investigation to two-component discs 
of stars and gas, and analysed the stability of THINGS spi- 
rals. They showed that ISM turbulence alters the condition 
for star-gas decoupling and increases the least stable wave- 
length, but hardly modifies the Q parameter at scales larger 
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Figure 4. The stability level of THINGS spirals. The diagnostics used in the left and right panels are our Q2 and Q3 parameters [see 
Eq. (19) for N = 2, 3]. The galaxy data are from Leroy et al. (2008), R is the galactocentric distance, and R25 is the optical radius. The 
data are colour-coded so as to show whether the stability level is dominated by stars or gas (H1/H2), as predicted by Eq. (19). In the 
right panel, the three data points that lie well below the critical stability level tell us that the nuclear region of NGC 6946 is subject 
to strong H2-dominated instabilities. This is consistent with the facts that NGC 6946 hosts a nuclear starburst (e.g., Engelbracht et al. 
1996) and a nuclear 'bar within bar' (e.g., Fathi et al. 2007). 
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Figure 5. Stability characteristics of THINGS spirals. The radial profiles shown in the left and right panels are computed by binning 
the Qjv data of Fig. 4 in 12 rings of width 0.1 R/R25, where Qjv is our Q stability parameter for multi-component and realistically thick 
discs [see Eq. (19) for = 2,3], i? is the galactocentric distance, and R25 is the optical radius. The right panel shows how frequently 
the stability level is dominated by stars or gas (H 1/H2 ) , and how frequently Q jv < 1 ■ The thick solid lines are the predictions of a 
three-component analysis (stars plus Hi plus H2), while the thin dashed lines represent the two-component case (stars plus gas). 



than about 100 pc. Since these are the usual scales of inter- 
est, we do not include that effect in our approximation. 

2.5 Application to THINGS spirals 

Our approximation is now complete. Let us then show how 
to use our Qn parameter for analysing the stability of galac- 
tic discs, and illustrate the strength of a multi-component 
analysis. 

In the following, we generalize the two-component ap- 



proach of Romeo & Wiegert (2011). We consider the same 
sample of spiral galaxies as in Sect. 2.1 and refer to Leroy 
et al. (2008), hereafter LOS, for a detailed description of the 
data and their translation into physical quantities. We treat 
stars. Hi and H2 as three components with the same surface 
densities and cr* as in LOS, but with distinct values of a^: 
(thi = 11 kms~^ (LOS), and aH2 = 6 kms~^ (Wilson et al. 
2011). For each galaxy, we compute the radial profile of Q3 
using Eq. (19). We adopt (cfz/crij)* = 0.6, as was assumed 
by LOS, and (ctz/ct j{)hi = {o'z/o-R)n2 ~ 1, as is natural for 



© 2013 RAS, MNRAS OOO.[HfT0l 



A simple and accurate approximation for Q 7 



coUisional components. For comparison purposes, we also 
compute the radial profile of Q2, treating the ISM as a sin- 
gle component with Eg — Ehi + Sh2 and erg = 11 kms~ 
(LOS). 

Fig. 4 shows Q3{R) and Q2{R) for the whole galaxy 
sample. Note that the Q3 data are characterized by a sharp 
transition at about half the optical radius. For 7i<0.6f?25, 
Q3 spans a range of one order of magnitude, and a sig- 
nificant fraction of the data lie below or near the critical 
stability level (although in most of the cases Q3 > 1). For 
R> 0.6i?25, Q3 varies within a narrow range of values, and 
there is a single data point with Qa < 1. Why do the inner 
and the outer discs of THINGS spirals have distinct stability 
properties? Why does the transition occur at about half the 
optical radius? To answer these questions, we have colour- 
coded the Q3 data so as to show which component domi- 
nates the local stability level. This is an important piece of 
information, which can easily be predicted using our Qjv di- 
agnostic [see Eq. (19)]. The fundamental difference between 
inner and outer spiral discs is how II2 contributes to disc 
(in)stability. For R25, H2 dominates in one third of 

the cases: it lowers the overall stability level, and increases 
the variance of Q3. At R ~ 0.6 R25, the contribution of 
II2 becomes negligible. Thus H2 leaves a characteristic im- 
print on the stability of THINGS spirals, even though stars 
dominate in most of the cases. The contribution of H I is in- 
stead negligible everywhere, even at the edge of the optical 
disc, where Hi is expected to contribute significantly. Such 
a stability scenario cannot be predicted by a two-component 
analyis. Note, in fact, that the Q2 data underestimate signif- 
icantly how gas contributes to disc (in)stability, and fail to 
reproduce the transition at half the optical radius (compare 
the left and right panels of Fig. 4). 

Fig. 5 illustrates how the stability properties of 
THINGS spirals vary with galactocentric distance. To ex- 
tract such information, we have binned the Q3 data of Fig. 
4 in 12 rings of width 0.1ii/i?25. For each ring, we have 
computed the median of Q3, the percentage of cases in 
which Qs < 1, and how frequently each component dom- 
inates the local stability level. For comparison purposes, we 
have also binned the Q2 data and computed the correspond- 
ing stability characteristics. Note that H2 plays a primary 
role in disc (in)stability for R < 0.1-0.2 R25, i.e. up to dis- 
tances of about one disc scale length (-R25 ~ 4.6 ± 0.8 Rd; 
LOS). Thereafter stars dominate more often. Note also that 
the frequency of H2-dominated cases decreases markedly 
with galactocentric distance, and falls below 10% a.t R = 
0.5-0.6 -R25. This corresponds to the transition radius iden- 
tified in Fig. 4, and to the outskirts of the expected H2 
domain (Eh2 > Em for R < 0.43 ± 0.18 R25; LOS). The fre- 
quency of H l-dominated cases shows a tendency to increase 
for R > 0.7-0.8 R25, but it never rises above 10%. The re- 
maining stability characteristics provide more intriguing in- 
formation. The frequency of locally unstable cases is below 
10%, except at distances smaller than about half the disc 
scale length. The median of Q3 lies well above the critical 
stability level, and is remarkably constant across the entire 
optical disc: QsmediR) — 2 (see the left panel of Fig. 5, and 
note the linear scale on the y-axis). Note, finally, how fast the 
two-component case diverges from the predicted radial pro- 
files as we approach the galactic centre (compare the dashed 
and solid lines in the two panels of Fig. 5). This result il- 
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Figure 6. Consistency between the results of our stability anal- 
ysis and the short- wavelength approximation, fcmaxfl 2> 1, where 
fcmax is the least stable radial wavenumber, R is the galactocen- 
tric distance, and R25 is the optical radius. Large error bars signal 
a transition from star- to H2-dominated regimes, or vice versa, 
which causes a jump in the value of fcmax (see Sect. 2). 



lustrates, once again, (i) how important it is to treat stars, 
H I and H2 as three distinct components when analysing the 
stability of galactic discs, and (ii) the strong advantage of 
using our Qjv parameter as a stability diagnostic. 



3 DISCUSSION 

Now that we have illustrated the strength of Qjv, let us re- 
member its weaknesses. Like all Q parameters and stability 
criteria discussed so far, Qn measures the stability of the 
disc against local axisymmetric perturbations, so it assumes 
that kR 3> 1. This is the short- wavelength approximation, 
which Binney & Tremaine (2008) define as 'an indispensable 
tool for understanding the properties of density waves in dif- 
ferentially rotating disks' (see p. 485 of Galactic Dynamics). 
Here the relevant k is the least stable radial wavenumber, 
which can be approximated as fcmax ~ i^/(^m (see Sect. 2). 

Rather than estimating the magnitude and radial be- 
haviour of kmnxR with qualitative arguments, we plot this 
quantity as a function of galactocentric distance for all 
THINGS spirals analysed in Sect. 2. Fig. 6 illustrates that 
the results of our stability analysis are consistent with the 
short- wavelength approximation. For instance, the condition 
fcmaxJ? > 3 is fulfilled in 93% of the cases, and is always true 
for i? > 0.25 R2B ~ ^d, i.e. at distances larger than about one 
disc scale length (remember that R25 = 4.6 ± 0.8 Rd; LOS). 
In contrast, there are only 2% of the data with kmaxR < 1, 
all close to the galactic centre: i? < 0.05 R25 ~ 0.25 Rd- Such 
data correspond to star-dominated regimes, which are more 
prone to violate the short-wavelength approximation since 
fcmax ~ k/ct* and a* > am > aH2- A comparison with 
the right panel of Fig. 4 shows that the corresponding val- 
ues of Qs are well above unity. As Q3 is well above unity 
and star-dominated in most of the cases, the few data with 
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k^BxR < 1 do not have a significant influence on tlie stability 
of THINGS spirals. Hi- and H2-dominated regimes are fully 
consistent with the short-wavelength approximation, and so 
are the corresponding values of Q3. In particular, the three 
data points that lie well below the critical stability level tell 
us that the nuclear region of NGC 6946 is subject to strong 
H2-dominated instabilities (see again the right panel of Fig. 
4). This is consistent with the facts that NGC 6946 hosts a 
nuclear starburst (e.g., Engelbracht et al. 1996) and a nu- 
clear 'bar within bar' (e.g., Fathi et al. 2007). 

While the short-wavelength approximation is satisfied 
by most spiral galaxies, the assumption of axisymmetric (or 
tightly wound) perturbations is not so general. Local non- 
axisymmetric stability criteria are far more complex than 
Toomre's criterion: they depend critically on how tightly 
wound the perturbations are, and cannot generally be ex- 
pressed in terms of a single 'effective' Q parameter (e.g., 
Lau & Bertin 1978; Morozov & Khoperskov 1986; Bertin et 
al. 1989b; Jog 1992; Lou & Fan 1998; Griv & Gedalin 2012). 
However, there is a general consensus that non-axisymmetric 
perturbations have a destabilizing efTect, i.e. a disc with 
Q > 1 can still be locally unstable against such perturba- 
tions. Gas dissipation has a similar efTect (Elmegreen 2011). 
These may be two of the reasons why the stability level of 
THINGS spirals is, on average, well above unity. The re- 
markable flatness of Q3mcd(J?) across the entire optical disc 
is far more intriguing. 

The assumption of local perturbations is quite contro- 
versial. While there is a general consensus that locally sta- 
ble discs can be globally unstable as regards spiral structure 
formation, the dynamics and evolution of spiral structure 
depend critically on the radial proflle of the Q stability pa- 
rameter (e.g., Bertin et al. 1989a, b; Lowe et al. 1994; Romeo 
1994; Korchagin et al. 2000, 2005; Khoperskov et al. 2007; 
Sellwood 2011; Khoperskov et al. 2012). Our results about 
the stability level of THINGS spirals have no direct impli- 
cations for that problem because they concern the THINGS 
sample as a whole, not each of the spirals. Useful constraints 
on the nature of spiral structure in galaxies might be found 
by analysing the radial profile of Qiv for each THINGS spi- 
ral, and by searching for trends in Qn{R) along the Hub- 
ble sequence. This is however well beyond the scope of the 
present paper. 



4 CONCLUSIONS 

This paper provides a simple analytical recipe for esti- 
mating the Q stability parameter in multi-component and 
realistically thick discs [see Eq. (19)]. Our Qjv parame- 
ter applies for any number of stellar and gaseous compo- 
nents (i = 1, . . . , A''), and for the whole range of velocity 
anisotropy observed in galactic discs: < {crz/crR,)i ^ 1. The 
accuracy of this approximation can be rigorously quantified 
in the thin-disc limit, where it scales as N~^^'^ . For N = 3, 
the predicted root-mean-square error is well below 10%. This 
is true even for larger values of A*'. For example, in the ten- 
component model (s) of the Solar neighbourhood considered 
by Rafikov (2001) our approximation is accurate to within 
5%. A further strength of the Qm diagnostic is that it pre- 
dicts which component dominates the local stability level. 



This is a useful piece of information, which should always 
be given when analysing the stability of galactic discs. 

This paper also provides the first three-component anal- 
ysis of THINGS spirals. Our analysis predicts how stars, H I 
and H2 contribute to disc (in)stability, and how the stability 
properties of such galaxies vary with galactocentric distance. 
We show that H2 plays a primary role up to distances of 
about one disc scale length. Stars dominate thereafter, but 
the contribution of H2 remains significant even at distances 
as large as half the optical radius. This is in sharp contrast 
to the role played by Hi, which is negligible up to the edge 
of the optical disc. We also show that the stability level 
of THINGS spirals is, on average, remarkably fiat and well 
above unity. 
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APPENDIX A: DERIVATION OF EQ. (17) 

From standard error analysis, we know that the relative un- 
certainty of Qn is approximately equal to that of 1/Qn- 



Table Bl. Rafikov's reference model of the Solar neighbourhood: 
surface densities and velocity dispersions of the various compo- 
nents. 



i 


Component 


L W y J 


fjj [km s 


1 


Tci\ /r 


13.0 


/ .0 


2 


giants 


0.4 


26.0 


3 


Mv < 2.5 


0.9 


17.0 


4 


2.5 <Mv < 3.0 


0.6 


20.0 


5 


3.0 <Mv < 4.0 


1.1 


22.5 


6 


4.0 <Mv < 5.0 


2.0 


26.0 


7 


5.0 <Mv < 8.0 


6.5 


30.5 


8 


Mv > 8.0 


12.3 


32.5 


9 


white dwarfs 


4.4 


32.5 


10 


brown dwarfs 


6.2 


32.5 



AQiv A(l/Qiv) 



(Al) 



Qjv 1/Qjv 

(see, e.g., Bevington & Robinson 2003). We also know that 
A(l/Qjv) arises from estimation uncertainties in Wi [see Eq. 
(15)]. If we assume that all AWi are uncorrelated and of 
comparable magnitude, then the error propagation equation 
reduces to 

/ N ^ \ 1/2 



A(l/Q^) ~ AW > ^ 



(A2) 



Remembering that Qi > Qm, we get 
A(l/Q. ' - 



< . 



(A3) 



Noting that 1/Qjv > 1/Qm and using Eq. (Al), we then 
find 

AQn 



<AW\^. 



(A4) 



AW can be evaluated from the two-component case, where 
the root-mean-square value of the relative error (Q2 — 
QRk)/QRk is 5% (see Fig. 2). Setting AW V2 ~ 0.05, we 
finally obtain AW ~ 0.03 and thus 



AQn 
Qn 



<0.03 VN 



(A5) 



APPENDIX B: THE TEN-COMPONENT CASES 
ANALYSED BY RAFIKOV (2001) 

• Case 1: Rafikov's reference model of the Solar neigh- 
bourhood. The epicyclic frequency is k = 36 kms~^ kpc~^. 
The surface densities and velocity dispersions of the various 
components are listed in Table Bl. 

• Case 2: same as Case 1, but with the velocity disper- 
sion of white and brown dwarfs decreased from 32.5 to 20.0 
kms"'^. 

• Case 3: same as Case 1, but with the total surface den- 
sity of white and brown dwarfs increased from 10.6 to 25.0 
Mq pc-2. 

• Case 4: same as Case 1, but with the gas velocity dis- 
persion decreased from 7.0 to 5.9 kms~^. 

• Case 4+1- same as Case 1, but with the gas surface 
density increased from 13.0 to 14.8 Mq pc~^ (this case was 
only mentioned by Rafikov). 
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